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The resilience to averaging over an initial energy distribu- 
tion of reducibility and thermal scaling observed in nuclear 
multifragmentation is studied. Poissonian reducibility and 
the associated thermal scaling of the mean are shown to be 
robust. Binomial reducibility and thermal scaling of the ele- 
mentary probability are robust under a broad range of condi- 
tions. The experimental data do not show any indication of 
deviation due to averaging. 
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The complexity of nuclear multifragmentation under- 
went a remarkable simplification when it was empirically 
observed that many aspects of this process were: a) "re- 
ducible" ; and b) "thermally scalable" |l|-^] . 

"Reducibility" means that a given many-fragment 
probability can be expressed in terms of a corresponding 
one- fragment probability, i.e., the fragments are emitted 
essentially independent of one another. 

"Thermal scaling" means that the one-fragment prob- 
ability so extracted has a thermal- like dependence, i.e., 
it is essentially a Boltzmann factor. 

Both "reducibility" and "thermal scaling" were ob- 
served in terms of a global variable, the transverse en- 
ergy E t (defined as E t — J^-Ejsin i.e. the sum of 
the kinetic energies E of all charged particles in an event 
weighted by the sine squared of their polar angles 9), 
which was assumed (see below) to be proportional to the 
excitation energy of the decaying source(s) 

In particular, it was found that the ^-integrated mul- 
tiplicity distributions P(n) were binomially distributed, 
and thus "reducible" to a one-fragment probability p. 
With higher resolution, it was noticed that for each in- 
dividual fragment species of a given Z, the n^-fragment 
multiplicities P(nz) obeyed a nearly Poisson distribu- 
tion, and were thus "reducible" to a single-fragment 
probability proportional to the mean value {nz) for each 

The one-fragment probabilities p showed "thermal 
scaling" by giving linear Arrhenius plots of hip vs \j\fE~t 
where it is assumed that \jE~ t cx T. Similarly n-fragment 
charge distributions P n {Z) were shown to be both "re- 
ducible" to a one-fragment Z distribution as well as 
"thermally scalable" H. Even the two-fragment angu- 
lar correlations Pi^A^) were shown to be expressible 
in terms of a one-body angular distribution with ampli- 
tudes that are "thermally scalable" Q. Table | gives a 
summary of the "reducible" and "thermal scaling" ob- 
servables. 
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TABLE I. Summary of reducible and thermal scaling ob- 
servables in nuclear multifragmentation. 



Empirically, "reducibility" and "thermal scaling" are 
pervasive features of nuclear multifragmentation. "Re- 
ducibility" proves nearly stochastic emission. "Thermal 
scaling" gives an indication of thermalization. 

Recently, there have been some questions on the signifi- 
cance (not the factuality) of "reducibility" and "thermal 
scaling" in the binomial decomposition of Z-integrated 
multiplicities Q. For instance, had the original distribu- 
tion in the true excitation-energy variable been binomi- 
ally distributed and thermally scalable, wouldn't the pro- 
cess of transforming from excitation energy E to trans- 
verse energy E t through an (assumedly) broad transfor- 
mation function P{E,Et) destroy both features? 

Specifically, under a special choice of averaging func- 
tion (Gaussian), for a special choice of parameters (vari- 
ance from GEMINI Q), and for special input p (the ex- 
citation energy dependent one-fragment emission proba- 
bility) and m (the number of "throws" or attempts) to 
the binomial function, the binomial parameters extracted 
from the averaged binomial distribution are catastrophi- 
cally altered, and the initial thermal scaling is spoiled M . 
This "spoiling" in is not due to detector acceptance 
effects (which has been commented on extensively in || ) , 
but rather is due to the intrinsic width of correlation be- 
tween Et and E as discussed below. 

It should be pointed out that, while the decomposition 
of the many-fragment emission probabilities P(n) into p 
and to may be sensitive to the averaging process, the 
quantity {mp) is not 0. However, both p and {mp) are 
known to give linear Arrhenius plots with essentially the 
same slope (see below) . This by itself demonstrates that 
no damaging average is occurring. 

Furthermore, we have observed that by restricting 
the definition of "fragment" to a single Z, the multi- 
plicity distributions become nearly Poissonian and thus 
are characterized by the average multiplicity {mp) which 
gives well behaved Arrhenius plots [Q . Thus, the linearity 
of the Arrhenius plots of both p and {mp) extracted from 
all fragments, and the linearity of the Arrhenius plots 
of {mp) for each individual Z value eliminate observa- 
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tionally the criticisms described above. In fact, it follows 
that no visible damage is inflicted by the true physical 
transformation from E to E t . Therefore, the experimen- 
tal Poisson "reducibility" of multiplicity distributions for 
each individual Z and the associated "thermal scaling" of 
the means eliminates observationally these criticisms. 

We proceed now to show in detail that: 1) binomial 
reducibility and thermal scaling are also quite robust un- 
der reasonable averaging conditions; 2) the data do not 
show any indication of pathological behavior. 

We first discuss the possible origin and widths of the 
averaging distribution. 

It is not apparent why the variance of P(E, E t ) calcu- 
lated from GEMINI f§ should be relevant. GEMINI is a 
low energy statistical code and is singularly unable to re- 
produce intermediate mass fragment (IMF:3 < Z < 20) 
multiplicities, the magnitudes of E t , and other multifrag- 
mentation features. There is no reason to expect that the 
variance in question is realistic. 

Apparently, E t does not originate in the late thermal 
phase of the reaction. Rather, it seems to be dominated 
by the initial stages of the collision. Consequently its 
magnitude may reflect the geometry of the reaction and 
the consequent energy deposition in terms of the number 
of primary nucleoli- nucleon collisions. This is attested to 
by the magnitude of E t which is several times larger than 
predicted by any thermal model. Thus, the worrisome 
"thermal widths" are presumably irrelevant. 

Since there is no reliable way to determine the actual 
resolution of the correlation between E t and E, exper- 
imentally or via simulation calculations Q|, instead of 
using large or small variances, we will show: 

a) which variables control the divergence assuming a 
Gaussian distribution, and in what range of values the 
averaging is "safe", i.e. it does not produce divergent be- 
havior; 

b) that the use of Gaussian tails is dangerous and im- 
proper unless one shows that the physics itself requires 
such tails. 

The input binomial distribution is characterized by m, 
the number of throws (assumed constant in the calcu- 
lations in j7|), and p which has a characteristic energy 
dependence of 
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We denote the one-to-one image of E in E t space with a 
prime symbol. 

The averaging in Q| is performed by integrating 
the product of an exponential folded with a Gaussian 
(Eq. (12) Of §). 
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If the slope of the exponential is large, there will be 1) a 
substantial shift e in the peak of the integrand, and 2) a 
great sensitivity to the tail of the Gaussian. 
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FIG. 1. Upper left panel: the distorted Arrhenius plot 
for a value of B=40MeV 1/ ' 2 , m=12, and a fixed ratio of 
r.E t /-Et=0.3. The open circles represent the apparent sin- 
gle fragment emission probabilities extracted after a folding 
of the "thermal" binomial emission probabilities with a Gaus- 
sian distribution (see Eqs. (j5|)-(^|)). The solid circles and open 
diamonds show the effect of truncating the Gaussian tails at 
2a and la, respectively. The star symbols demonstrate results 
from folding with a square distribution of the same full width 
at half maximum as the Gaussian. The solid line represents 
the undistorted Arrhenius plot. The dashed line represents 
an experimental "soft" limit beyond which the transverse en- 
ergy may be unreliable as a measure of impact parameter |J 
or deposited energy. Lower left panel: same as upper left but 
for r E t I Et=Q.2. Upper right panel: the distorted values of 
the number of throws m apP as a function of E t for a value of 
B=40MeV 1/2 and a fixed ratio of r Et /E t =Q.3. The solid line 
represents the undistorted value of m=12. Lower right panel: 
same as upper right but for r E t / Et— 0.2. 
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This illustrates the divergence at small values of xq both 
in the shift of the integrand in (p) and (p 2 ) and the 

corresponding divergence in a 2 = (p 2 ) — (p) 2 . The scale 
of the divergence is set by the product a 2 B. Thus one 
can force a catastrophic blowup by choosing a large value 
of a 2 , of B, or of both. This is what has been shown to 
happen with large values of a 2 and B. The counterpart 
to this is that there possibly exists a range of values for 
B and a 2 which leads to a "safe" averaging process. 
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In order to illustrate this, we have calculated the "ap- 
parent" values of the single fragment emission probabil- 
ity Papp for widths characterized by the ratio of the full 
width at half maximum T E t over E t . Specifically we have 
extracted 



Papp 
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and m n 



L app 



Papp 

by calculating the observed mean: 



and variance: 
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for "thermal" emission probabilities P™ folded with a 
Gaussian distribution g(E t ). We have assumed to is con- 
stant. 

For a value of T Et /E t ^0.3, m=12, and P^TOMeV 1 / 2 
(consistent with the upper limits of the slopes observed in 
the Xe induced reactions @|§), the onset of divergence is 
observed in the Arrhenius plot at small values of E t (top 
left panel of Fig. |l], open circles). For r 'g t / 'E t —0.2 (open 
circles in bottom left panel of Fig. |l]), the divergent be- 
havior is "shifted" to even lower energies and the result- 
ing Arrhenius plot remains approximately linear. There- 
fore, the thermal signature survives. For both widths, 
the linear (thermal) scaling survives in the physically ex- 
plored range of l/\/E~t < 0.08 (E t > 150 MeV) shown 
by the dashed lines in Fig. [|. As we shall see below, 
the effect is weaker for even lower values of B which are 
commonly seen experimentally. 

The divergent behavior manifests itself as well in the 
parameter to, the number of "throws" in the binomial 
description. Values of m app are plotted (open circles) 
as a function of E t in the right column of Fig. |l| for 
T E J Et—O.S (top panel) and T Et /E t =0.2 (bottom panel). 

While the distortions depend mostly on the variance of 
the energy distribution, distributions with similar widths 
can be associated with very different variances. For in- 
stance, a Lorentzian distribution with finite T has infinite 
variance. Its use would lead to a divergence even for in- 
finitely small values of T. Thus, even innocent trimmings 
to the (non-physical) tails of a Gaussian can produce big 
differences in the variance of the distribution and in the 
ensuing corrections. We exemplify this point in two ways. 

a) We use a "square" distribution with a width equal 
to the full width at half maximum of the Gaussian. As 
can be seen by the star symbols of Fig. |l| this simple 
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FIG. 2. Top panel: the divergence energy (the energy at 
which m ap p and p apP change sign) as a function of the slope 
parameter B for a fixed ratio of T E t / Et=0.?>. Bottom panel: 
the divergence energy as a function of the width TE t /Et for 
a fixed value of B=40MeV 1/2 . 



exercise dramatically extends the range over which the 
average can be performed safely. 

b) We truncate the tails of the Gaussian at la (dia- 
monds) and 2<T (solid circles) in Fig. [I]. Already the cut 
at 2(7 shows a dramatic improvement over a full Gaus- 
sian. The la cut actually makes things even better than 
the square distribution (as seen in Fig. Q). 

To illustrate the conditions under which the "thermal" 
scaling survives (i.e. linear Arrhenius plots as a func- 
tion of 1/ \fE~t) , we have traced the evolution of the "di- 
vergence energy" (or the point at which m avv and p app 
change sign) as a function of the two parameters which 
control the strength of the divergence: the slope parame- 
ter B and the variance a 2 (hereafter characterized by its 
full width at half maximum value r# t « 2.35'\/o^). 

A particular example for r^ t /£ , t =0.3 is shown by the 
open circles in the top panel of Fig. ^. In addition, values 
of the divergence energy for la and 2a truncations of the 
Gaussian as well as a square distribution are also plotted. 
For all intents and purposes, divergencies that occur at 
less than 100 MeV do not alter substantially the linear 
Arrhenius plots as they have been observed to date 
in the E t range of 150 to 1600 MeV. 

In a similar manner, the dependence of the divergence 
energy can also be determined as a function of the relative 
width rE t /Et (for a fixed value of B). This behavior is 
demonstrated in the bottom panel of Fig. ||. 

A more global view of the parameter space is shown 
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FIG. 3. The divergence energies (contour line) as a function 
of the width TE t /Et and the slope B for m=12 and a Gaussian 
distribution truncated at 2a. 

in Fig. H where the divergence energy is plotted (contour 
lines) as a function of the width TE t /Et and slope B. The 
shape of the contour lines reflects the a 2 B scale deduced 
in Eqs. (^) and (|). The calculation in ref. [?J sits nearly 
in the upper right hand corner of the graph. But, as is 
clearly demonstrated, large regions exist where binomial 
rcducibility and thermal scaling survive (roughly given by 
the region with divergence energies less than 100 MeV). 

From the above exercises it is concluded that there is 
abundant room for the survival of binomiality and ther- 
mal scaling. 

In this second part, we show that none of the symp- 
toms of divergence are present in the available experi- 
mental data jl]-^. Furthermore, the average fragment 
multiplicity (n) is expected to be "distortion free" |7j. 
As such, it provides a baseline reference with which to 
compare the "distorted" variable, p app (to verify whether 
the label "distorted" is appropriate). In addition, we can 
force the divergence to appear in the data, by artificially 
broadening the E t bins, thus establishing that it is not 
present with ordinary (small) E t bins. Finally, we show 
that thermal scaling is present and persists in the data 
even when the divergence is forced. 

First, we draw attention again to the two pathologic 
features arising from excessive averaging. 1) The quan- 
tity m diverges near E t =0. 2) The quantity 1/p suffers 
a corresponding discontinuity at the same low energy. 

Inspection of the published data shows that: 

1) m never diverges near E t = 0. To the contrary 
m remains relatively constant or actually decreases with 
decreasing E t . This is particularly true for all of the Xe 
induced reactions [||,|| (see Fig. g); 

2) log 1/p is nearly linear vs. \/^/El over the experi- 
mental E t range without the indications of trouble sug- 
gested by the calculations in the previous section. 

Thus the experimental data do not show any signs of 
pathological features. 

The quantity (n) = m app p app does not suffer from the 
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FIG. 4. The inverse of the single fragment emission prob- 
ability (solid circles) and the inverse of the average fragment 
multiplicity (open circles) as a function of l/\/E~l for the re- 
action Ar+Au at E/A=11Q MeV. The solid lines are linear 
fits to the data. 



distortions due to averaging. In fact, (n) is a suitable 
alternative for constructing an Arrhenius plot in those 
cases where m depends only weakly on E t (as observed 
in many of the data sets we have studied) . A comparison 
of the Arrhenius plots constructed from 1 /p app and 1 / (n) 
is shown in Fig. m The striking feature of this compari- 
son is that the l/p apP values have the same slope as the 
"distortion-free" case of 1/ (n) . Similar observations can 
be made for all the other reactions studied so far. As a 
consequence both the "fragile" p and the "robust" (mp) 
survive the physical transformation P(E,E t ) unscathed. 

When the probability becomes small, the binomial dis- 
tribution reduces to a Poisson distribution. This can be 
achieved experimentally by limiting the selection to a sin- 
gle Z 0| . The observed average multiplicity is now exper- 
imentally equal to the variance. Thus we are in the Pois- 
son reducibility regime and can check the thermal scal- 
ing directly on (nz). For a Poisson distribution, log (nz) 
should scale linearly with \j\fE~t. This can be seen ex- 
perimentally for the average yield of individual elements 
of a given charge (see Fig. ||) for the reaction Ar+Au at 
_E/j4=110 MeV. For the case of a single species, the re- 
ducibility is Poissonian, and the thermal (linear) scaling 
with \j\fE~t is readily apparent. As pointed out at the 
outset of the paper, this evidence, together with that of 
Fig. ^ indicates that no significant averaging is occurring 
even in the case of binomial decomposition. 

The data can be "encouraged" to demonstrate the sort 
of catastrophic failures described here. By widening the 
bins in transverse energy (AE t ), we can induce an arti- 
ficial broadening to mimic a broad correlation between 
E and E t . For example, the behavior of p app and m app 
is shown in Fig. ^ for three different widths and two dif- 



ferent reactions. The divergencies of p app and m app are 
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FIG. 5. The average yield per event of elements with Z 
between 3 and 12 as a function of 1/s/Et. The lines are linear 
fits to the data. 
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readily visible for large AE t values, but are noticeably 
absent for small values. The spectacularly large bin- 
ning in E t (100 MeV!) necessary to force the anticipated 
pathologies to appear is reassuring indeed. Notice that 
here the absolute width, not the relative width, was kept 
fixed even at the lowest energies! Furthermore, the sta- 
bility of (n) is readily apparent from the complete overlap 
of the values of (n) extracted for different windows of E t 
(open symbols of Fig. ^ . 
In summary: 

a) Binomial reducibility and the associated thermal 
scaling survive in a broad range of parameter space. The 
single case shown in Q is an extreme one based on un- 
supported assumptions about the averaging function. 

b) The experimentally observed simultaneous survival 
of the linear Arrhenius plot for parameter p and the ro- 
bust average (mp) suggests that no serious damage is 
generated by the physical transformation P(E, E t ). 

c) The multiplicity distributions for any given Z value 
are Poissonian and the resulting average multiplicity 
(n) = (mp) gives linear Arrhenius plots confirming the 
conclusion in b). 

d) Finally, the data themselves do not show any in- 
dication of pathological behavior. This can be seen, for 
instance, by comparing the behavior of p with (n) . The 
pathology can be forced upon the data by excessively 
widening the E t bins. Even then, the thermal scaling 
survives in the average multiplicity. 
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FIG. 6. Left panels: the inverse of the single fragment emis- 
sion probability (solid symbols) and the inverse average frag- 
ment multiplicity (open symbols) as a function of 1/y/Et for 
the indicated bin widths AE t and systems. Right panels: the 
observed values of m as a function of Et for the indicated bin 
widths AEt and systems 
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